Electromagnetically Induced Transparency with Rydberg Atoms 



O 

(N 
-i— > 
O 

O 

OS 



X3 
Oh- 

> ! 

c : 

=5 

cr 



> 
O 
SO 

m 

SO 

o 



X 



David Pctrosyan, 1,2 Johannes Otterbach, 1,3 and Michael Fleischhauer 1 

1 Fachbereich Physik und Forschungszentrum OPTIMAS, 
Technische Universitat Kaiserslautern, D-67663 Kaiserslautern, Germany 
2 Institute of Electronic Structure and Laser, FORTH, GR-71110 Heraklion, Crete, Greece 
3 Physics Department, Harvard University, Cambridge, Massachisets 02138, USA 

We present a theory of electromagnetically induced transparency in a cold ensemble of strongly 
interacting Rydberg atoms. Long-range interactions between the atoms constrain the medium to 
behave as a collection of superatoms, each comprising a blockade volume that can accommodate 
at most one Rydberg excitation. The propagation of a probe field is affected by its two-photon 
correlations within the blockade distance, which are strongly damped due to low saturation threshold 
of the superatoms. Our model is computationally very efficient and is in quantitative agreement 
with the results of recent experiment of Pritchard et al. [Phys. Rev. Lett. 105, 193603 (2010)] 

PACS numbers: 42.50.Gy, 42.65.-k, 32.80.Ee, 



Strong dipole-dipole or van der Waals (VdW) inter- 
actions between atoms in highly excited Rydberg states 
[l| constitute the basis for promising quantum informa- 
tion schemes 0] and interesting many-body effects 0- 
HJ. Many of these studies utilize the dipole blockade 
mechanism which suppresses multiple Rydberg ex- 

citations within a certain interaction (blockade) volume. 
Electromagnetically induced transparency (EIT) [l4| can 
translate the interactions between Rydber g at oms into 
sizable interactions between single photons |15l4l7| . 



Recently, several experiments on EIT [18H21| , and the 
closely related CPT (coherent population trapping) [22| , 
with Rydberg atoms were performed. Strong VdW inter- 
actions between the atomic Rydberg states were promi- 
nently manifest in Ref. [21j: Increasing the probe field 



amplitude led to reduction of its transmission within the 
EIT window, which, quite surprisingly, was accompa- 
nied by negligible broadening and indiscernible shift of 
the EIT line. Here we develop a theoretical model for 
EIT with Rydberg atoms, whose predictions fully re- 



produce the experimental observations [211 ]. The crux 
of our approach is the coarse-grained treatment of the 
atomic medium composed of effective superatoms (SAs), 
with each SA represented by collective states of atoms 
in the blockade volume that can accommodate only one 
Rydberg excitation. A weak probe field propagates 
through the EIT medium with little attenuation, but for 
a stronger field with more than one photon per SA, the 
excess photons are subject to enhanced — essentially two- 
level atom — absorption. This leads to the field attenua- 
tion with the simultaneous buildup of an avoided volume 
between the probe photons 12 1. The inclusion of two- 



photon correlations is the key feature of our work, not 
present in the numerical simulations of [21( and recent 
theoretical studies Q which agreed with the experiment 
at weak probe fields but had significant discrepancies for 
stronger fields. Our theory is not limited to weak probe 
fields and/or low atomic densities, yet, despite intrin- 
sic nonlinearity, it is intuitive and numerically efficient, 



amounting to the solution of a pair of coupled differen- 
tial equations for the probe field intensity and its second- 
order correlation, in the spirit of the BBGKY hierarchy. 

Consider an ensemble of N = J v d 3 r p(r) cold atoms 
of density p(r) in the (quantization) volume V interact- 
ing with two optical fields. The quantized probe field £ p 
of frequency lo p acts on the atomic transition between 
the ground \g) and excited |e) states, and the control 
field of frequency ui c drives the transition |e) — > \r) with 
Rabi frequency f2 c [see Fig. HJa)]. A pair of atoms i 
and j at positions r< and rj excited to the Rydberg 
states \r) interact with each other via a VdW poten- 
tial hA(ri — Tj) = hC 6 /\ri — r\,| 6 [24J . In the frame 
rotating with frequencies w PjC , the system Hamiltonian 
% = %& + V a f + VvdW contains the unperturbed atomic 



part, % & = -h Y% [\al e + (A p - 
field and VdW interactions, V a f 



8 c )d J rr ], and the atom- 



R(l>£(2>> 



(a) 



■Vvaw 



10 '* < i « 

~i""f " c 



(b) r,([. r (l. 



IR (1, > 



\g>- 



Lis 



t p 




l£ (3) > 



-2)3£l p 



l£ UJ > 



IG> 



FIG. 1. (a) Level scheme of atoms interacting with probe Q p 
and control Sl c fields with detunings A p and 5 C . T e and r r 
are (population) decay rates of states \e) and \r), and VvdW 
denotes VdW interaction between atoms in Rydberg state 
\r). (b) Truncated level scheme of superatom, composed of 
n sa atoms, with corresponding transition amplitudes due the 
probe and control fields. 
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n c CT^ e 4-H.c] and Vvdw = hJ2i<j ^Afr.-rj)^, where where £^(17) is the projector onto the Rydberg exci- 
&i v = \fi)jj{v\ is the transition operator for atom j at tation of SA at 17. The physical meaning of the first 

and 5 C = lu c — u re are the term on the rhs of Eq. ((2]) is that an excited SA at 17 ~ r 

[T, RR (r) — > 1] induces divergent VdW shift averaged over 
the SA volume: A ~ ^- J y A(r')<fV -> oo. Actually, 
for a small cut-off in the interatomic separation, |A| ^> j e 
is finite but very large, which is the only relevant prop- 
erty. The last term s(r) = A(r - 17) £^(17) de- 
scribes the VdW shift induced by the external SAs out- 

(r) 

side the volume Vsa centered at r. It can be evaluated 
by replacing the summation by an integration over the 
entire volume V, excluding the SA at r, which, upon 
using the mean-field approximation, yields a small shift 
(S(r)) = f(t RR (r)). 



position Tj, A p = uj p — 
detunings of the probe and control fields, and fl p = rj£ p 
is the operator of the probe Rabi frequency, with r\ = 
p ge ^Ju) p / {2Tie$V) the atom- field coupling strength on the 
\g) — > \e) transition with dipole moment p ge . 

We consider the evolution of the probe field only 
along its propagation z axis, and assume uniform, un- 
depleted control field O c . Using Hamiltonian H, we ob- 
tain Heisenberg-Langevin equations for the field £ p {r) 
and continuous, appropriately averaged, atomic ^^(r) 
operators: 



(d t + cd z )£ p {r) = ir)Na ge (r), (la) 
d t a ge (r) = (iA p - j e )a ge (r) + itt* c a gr (r) 

+ih p (r)[a gg (r) - a ee (r)], (lb) 
d t a gr (r) = [i(A 2 - S(r)) - 7 r ]<T sr (r) 

-i£lp(r)a er {r) + ifl c a ge (r), (lc) 

where j e > ^T e and 7 r (<C 7 e ) are the transversal relax- 
ation rates, with the associated noise operators dropped, 
A 2 = A p + S c is the two-photon detuning, and S(r) = 
J dVp(r') A(r-r')oy r (r') is the total VdW induced shift 
of level \r) for an atom at position r. Since S(r) involves 
integration over all spatial coordinates r' £ V, Eqs. ([1]) 
are highly nonlocal. We therefore need to contrive an 
efficient method to evaluate the VdW shift S(r). 

We shall be concerned with stationary interaction and 
drop in Eqs. ([I} all the time-derivatives. Consider for 
the moment Eqs. (fTb|) and (|Tc|) without the VdW shift S 



and small relaxation 7 r terms. When O p . c , |A Pj 2| < 7e, 
we can approximate the population of Rydberg state \r) 
by a Lorcntzian function of A2: (6y r ) ~ (OIOp) /(|O c | 2 + 
A|7^/|^ c | 2 ), with the half- width w = \£l c \ /je- Observe 
now that an atom in Rydberg state \r) would induce 
VdW shift A(i?) of level \r) for another atom separated 
by distance R, which effectively translates into the two- 
photon detuning A2. The VdW interaction then blocks 



the excitation of all the atoms for which A(i?) 



> 



w. 



This 



is the essence of the Rydberg blockade 
define the blockade radius i? s 
ensemble of n sa = pV st 



We therefore 



^/Ce/w and call the 
atoms within volume V Ba = -j-i?!? a 
"superatom" (SA). Since in general the atomic density 
p(r) varies with position r, so does n sa ( r), but the density 
of SAs p sa = V^ 1 = ^ v /|rj c | 2 /(7 e C 6 ) is constant. 

Each SA can contain only one Rydberg excitation de- 
localized over V sa . We may therefore treat the medium 
as a collection of N sa = p sa V SAs at positions 17 , which 
implies a spatial coarse-graining with the grain size 2i? sa 



25|,|26|. The total VdW shift 5(r) at position r can then 



be expressed as 

iV sa 

£(r)«£A( 

j 



T^tmiiTj) = A± RR (r) + s(r), (2) 



Assuming that the probe field O p varies little over 
distance ~ R sa , we can describe the dynamics of in- 
dividual SAs in terms of collective states and oper- 
ators defined within the blockade volume V sa . The 
level scheme of SA is shown in Fig. QJb): \G) = 
\gi, g2, ■ ■ ■ , 9n Ba ) is the ground state, and \R^) = 



_1 yir 



151,52, 



■ 9r> 



is the single collective 



Rydberg excitation state, while are the prop- 

erly symmetrized (Dicke) states with k atoms in |e). 
The corresponding transition amplitudes {E^> | V a f \G) = 
y/T^Clp, V af = Cl c , etc., depend of the num- 

ber of atoms n aa in V sa . In order to calculate T, RRl we 
now proceed along the lines similar to the single atom 
treatment. Starting with the SA in \G), we adiabati- 
cally eliminate all the excited states \E^) having large 
widths ~ kje- Note that state |.R^) is reached from 
\G) by 2-photon transition, while all the other states 
require 2+k photon transitions; therefore their 
adiabatic elimination affects little the populations of \G) 
and \R^). We then obtain for the SA operators T<gr = 
\G)(R^\= n c ^Zh p t GG /[{Ap + i le )A 2 - \tt c \ 2 ] and 
T, RR = Yi R qYiq R . To account for possible saturation 
of transition \G) — > \R^) when the number density of 
probe photons p p hot is comparable to, or larger than, the 
density of SAs p sa , we take Egg + Err = 1, which finally 
yields [27| 



10, 



■•RR 



2 n Qt6 
ft sa J L p^ L p 



I O c \ 2 n sa ClpClp 



|O c P-A p A 2 p + Al7| 



(3) 



We next examine the probe field propagation in the 
atomic medium. For moderate Rabi frequency fl p < "f e 
and number density of photons p p hot *C p, we can assume 
linear response of individual atoms to the applied field, 
setting CT ee , o'er - > and a gg = 1. We then arrive at 
the propagation equation for the probe field amplitude, 
d z £ P = i^&Sp, where k = soP is the resonant (intensity) 
absorption coefficient proportional to the atomic absorp- 
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tion cross-section = w p |p ge | 2 / (heocj e ) , while 



&(r) = S fl i?(r) 



He 



le - «A p 



[1 - Sflfl(r)] 



'lie 



le - iA p + |fi c | 2 [ 7r - i(A 2 - (s(r)))] 



T (4) 



is the operator-valued polarizability. Here the first frac- 
tion is the polarizability axLA of a two-level atom, while 
the second fraction, barring the small shift (s(r)), is the 
usual EIT polarizability «eit [3] ■ Physically, if the SA 
at position r contains a Rydberg excitation [S^^(r) — > 
1], the two-photon detuning is shifted out of the EIT 
transparency window [A 3> j e ] and the probe field £ p 
sees an absorbing two-level system; if no Rydberg excita- 
tion is present, the medium response is that of usual EIT 
with a small mean-field shift due to the VdW interaction 
with the external SAs. Then the expectation value of the 
probe field intensity obeys the equation 

z (£t(r)4(r)) = - K (T)(8l(r)Im[a(r)}E p (v)). (5) 

Note that factorizing out Im[(d(r))] in a mean-field sense 
would amount to ne glec ting the essential two-particle 
quantum correlations 17j originating from nonlinear re- 
sponse of the atoms to the Rydberg excitations. We 
therefore proceed more carefully and replace &(r) in 
Eq. ([5]) by its expectation value conditioned upon the 
presence of photon at r, denoted by (-) r , 

(&(r)) r = (E Kii (r)) r aTLA + [1 - (t RR (r)) r }a E iT- (6) 

The conditional Rydberg population (S_r_r),. of the 
SA at r is obtained from Eq. ([3]) by the replace- 
ment {ll(r)n p (r) -> (f2t( r )n p (r)) g ( P \v), where the 

(2) 

probe field intensity correlation function g p (r) = 

di^iSISw) quantifics thc probability of hav - 

ing simultaneously at least two photons in the blockade 

(r) 

volume Vsk . Thc field intensity is now coupled to its 
two-photon correlation g p (r) which in turn evolves upon 
propagation. Note that linear, e.g. bare EIT, response 
of the medium does not change the correlation function 
of the propagating field, and only nonlinear, i.e. con- 

(2) 

ditional, absorption cx Im[(a(r)) — ckeit] modifies g p , 
which therefore obeys the equation of motion [28| 

d z9p 2 \v) = -K(r)<E JUl (r))Im[oTLA - a ElT ]g p 2 \r). (7) 

Hence, within the EIT transparency window, where 
Im[aEiT] — 0, the correlations between the photon pairs 
with relative distance smaller than the blockade (SA) ra- 
dius decay with the rate proportional to the probabil- 
ity of SA excitation (Err) and the absorption rate of a 
two-level system Im[o?TLA] ■ We note that our treatment 
involves only single transverse mode of the probe field, 
which is effectively defined by the SA cross-section. If, 
however, during propagation there is strong mixing of 




A/231 (MHz) 

FIG. 2. Top: Probe field transmission I p (L)/I p (0) versus 
detuning A p , for different input intensities corresponding to 
Q p (0)/2n = 0.15,0.5, 1.0 MHz. Thin lines are experimental 
curves from [2l[ , thicker lines are obtained via stochastic sim- 
ulations of Eqs. (J5j)-(JTI) averaged over 10 independent real- 
izations. Bottom: The corresponding intensity correlation 
functions g p (L). 



the transverse modes, it would preclude the buildup of 
(anti)corrclations between the photons. 

Given the input field "intensity" I p = (Cl p tl p ) and 

(2) 

its correlation function g p ' [for "classical" coherent field 

(2) 

9j> (0) = 1J) we then use thc following stochastic pro- 
cedure to spatially integrate the coupled coarse-grained 
Eqs. ©-0 for z <G [0,L]: We divide the propaga- 
tion distance L into L/(2R sa ) intervals corresponding to 
SAs, and for z within each SA we determine via Monte- 
Carlo sampling of (Ejjfl(r)) r whether the SA is excited, 
S_Ri?(r) — > 1, or not, Tjrr(y) 0. We then average over 
several independent realizations. The limit of infinitely 
many such realizations corresponds to continuous polar- 
izability of Eq. ©. 

We employ our theory to simulate the experiment of 
Ref. [2l[ with an ensemble of cold 87 Rb atoms: \g) = 
5S 1/2 \F = 2,m F = 2), |e) = 5P 3/2 \F = 3,m F = 
3) with T e = 3.8 x 10 7 s" 1 , and |r) = 60S 1/2 with 
IV = 5 x 10 3 s" 1 and C 6 /2tt = 1.4 x 10 11 s-^im 6 
[24j corresponding to repulsive VdW interactions. ^ e%r 
also include the one- and two-photon laser linewidths 
Slji^/Zk — (5.7,11) x 10 4 s _1 . The atomic density is 
p(z) = po ex p[— — z o) 2 /2<Tp] with peak po = 1.32 x 
10 7 mm~ 3 and half- width a p = 0.7mm; indistinguishable 
results are obtained for homogeneous ensemble of density 
p = 1.2 x 10 7 mm -3 and length L = 1.3 mm, leading to 
the resonant optical depth of kL = 4.524. The control 
field n c /2n = 2.25 x 10 6 s^ 1 ^ is slightly detuned by 
S c /2n = — 10 5 s -1 . The corresponding blockade radius is 
i? sa ~ 6.6pm and each SA contains on average n sa ~ 14.7 
atoms. We emphasize that our simulations are insensi- 
tive to moderate variations (±20%) of the SA volume V^ a 
and thc number of atoms n sa (~ 14 ± 3) it contains. 

In Fig. [5] we compare the transmission spectra for 
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FIG. 3. (a) Transmission (top) and intensity correlation (bot- 
tom) spectra of the probe field, for various input intensities, 
(b) The same, but setting gf\z) = lVz G [0,L]. (c) Peak 
probe transmission around the EIT line center, (d) EIT 
linewidth <5u;eit (FWHM) , and (e) detuning A™ ax at the max- 
imum T max , versus the input probe Rabi frequency fi p (0). 
The black lines (circles) correspond to case (a) and the brown 
lines (squares) to (b). A p ,5ojeit, 



different input probe intensities with the corresponding 
plots of Ref. dlj]. Already for ft p (0)/27r > 0.1 MHz the 
VdW interaction induced nonlincarities play an impor- 
tant role. The agreement between our stochastic simu- 
lations and the experiment is remarkable. We also show 

(2) 

the local intensity correlation g p ' (L) at the exit from the 
medium. 

Figure [3] summarizes the results of our simulations in- 
volving the continuous polarizability of Eq. ((6]). The 
weak field of O p /2-7r < 0.01 MHz encounters linear EIT 
response. Increasing the input probe intensity leads to 
lesser transmission through the EIT window (A2 ~ 0) 
and to small mean-field shift and broadening of the EIT 
line. This is due to the higher probability of two or 
more photons, exciting Rydberg states |r), to be at 
the same SA. The induced large VdW level shift A re- 
sults in strong photon absorption, simultaneously reduc- 
ing the photon coincidence probability within the SA 

(2) 

volume K,a- Hence, both I p (z) and g p (z) decay, but 
1, the attenuation of the probe field in- 



(2) 

once gp (z) -C 



v = 2\Q C \ 2 / (nj e ) (~ 6000 m/s) is the probe group ve- 
locity within the EIT window IA2I ^ Sloeit, we have 
that pphot = (p/4)(^p^p)/|^ c | 2 and the maximal satu- 
ration intensity is 7 p nax ~ (Ap sa / p)\Vl c \ 2 . In the medium 
the photons are anticorrelated (antibunched) within the 
temporal window of St ~ 2R aa /v (~ 1.6 ns), which does 
not change when they leave the medium for free space. 

Had we not taken into account the probe field intensity 

(2) 

correlation, equivalent to setting g p ' (z) = IV z € [0, L], 
Fig. [3jb), we would have faster, exponential decay of 
I p {z), unrestrained by the buildup of avoided volume be- 
tween the photons, as well as sizable shift and broadening 
of the EIT line, which contradict the observations of 21 



tensity I p (z) slows down. Eventually I p = (DlClp) sat- 
urates at a value corresponding to less than one pho- 
ton per SA, /Ophot Psa, with vanishing coincidence 
probability. With p p hot = ^ocI p / {2p 2 ge ui p v ), where 



Outside the EIT window, around the Autler-Towncs 
doublet A2 ~ ±fi c , the probe is strongly absorbed, 
Im[(a)] ~ 1, but the correlation function is amplified, 
since in Eq. ([7]) Im[cvxLA — «eit] < 0. In other words, 
linear absorption is larger than the conditional absorp- 
tion, which results in photon bunching but very low flux. 

We finally note that for relatively strong input fields 
^p(0) ^ 7e of Fig. [3l the validity of linear response of the 
atoms inherent in polarizability of Eq. (j4]) may not be a 
priori justified. This is indeed the case for an optically 
thin atomic medium. But in the optically thick medium, 
within a few absorption lengths, even a strong probe field 
and its two-photon correlation function quickly decay to 
the level at which the above approximation is justified. 

To conclude, EIT via atomic Rydberg states is sup- 
pressed by collective Rydberg excitations of SAs which 
depend on the local probe field intensity and its two- 
particle correlation within the SA (blockade) volume. 
For strong input fields, the buildup of anticorrclations 
between the photons upon propagation through the 
medium leads to the saturation of transmitted field inten- 
sity to a value corresponding to one photon per blockade 
volume. Conversely, suitably antibunched input fields 
should exhibit large transmission affected only by small 
linear absorption. 

In a one-dimensional configuration, the spatial corre- 
lations between the photons in the medium translate at 
the output into temporal correlations in free space, which 
can be measured by coincidence detection. The limit 
of maximal saturation intensity /™ ax of the transmitted 
through the EIT window field then corresponds to a train 
of non-overlapping single-photon pulses with the tempo- 
ral separation St of a few ns. 
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Gorshkov, Th. Pohl and M.D. Lukin for stimulating dis- 
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Foundation (D.P.), the Harvard Quantum Optics Center 
(J.O.), and SFB TR49 (M.F.). 



5 



[1] T.F. Gallagher, Rydberg Atoms (Cambridge University 

Press, Cambridge, 1994). 
[2] M. Saffman, T.G. Walker, and K. M0lmer, Rev. Mod. 

Phys. 82, 2313 (2010). 
[3] H. Weimer et al, Phys. Rev. Lett. 101, 250601 (2008). 
[4] J. Honer et al, Phys. Rev. Lett. 105, 160404 (2010). 
[5] H. Weimer and H.P. Biichler, Phys. Rev. Lett. 105, 

230403 (2010). 

[6] J. Schachenmayer et al, New J. Phys. 12, 103044 (2010). 
[7] T. Pohl, E. Dernier, and M.D. Lukin, Phys. Rev. Lett. 

104, 043002 (2010). 
[8] M. Viteau et at, Phys. Rev. Lett. 107, 060402 (2011). 
[9] M.D. Lukin et al, Phys. Rev. Lett. 87, 037901 (2001). 
[10] T. Vogt et al, Phys. Rev. Lett. 97, 083003 (2006). 
[11] D. Tong et al, Phys. Rev. Lett. 93, 063001 (2004). 
[12] K. Singer et al, Phys. Rev. Lett. 93, 163001 (2004). 
[13] R. Heidemann et al, Phys. Rev. Lett. 99, 163601 (2007). 
[14] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. 

Mod. Phys. 77, 633 (2005). 
[15] I. Friedler et al, Phys. Rev. A 72, 043803 (2005); B. He 
et al, ibid 83, 022312 (2011); E. Shahmoon et al, ibid 
83, 033806 (2011). 
[16] D. Petrosyan and M. Fleischhauer, Phys. Rev. Lett. 100, 
170501 (2008). 

[17] A.V. Gorshkov et al, Phys. Rev. Lett. 107, 133602 



(2011). 

[18] A.K. Mohapatra, T.R. Jackson, and C.S. Adams, Phys. 

Rev. Lett. 98, 113003 (2007). 
[19] A.K. Mohapatra et al, Nature Phys. 4, 890 (2008). 
[20] A. Tauschinsky et al, Phys. Rev. A 81, 063411 (2010). 
[21] J.D. Pritchard et al, Phys. Rev. Lett. 105, 193603 

(2010) . 

[22] H. Schempp et al, Phys. Rev. Lett. 104, 173602 (2010). 

[23] C. Ates, S. Sevingli, and T. Pohl, Phys. Rev. A 83, 
041802(R) (2011); S. Sevingli et al, Phys. Rev. Lett. 107, 
153001 (2011); S. Sevingli et al, J. Phys. B 44, 184018 

(2011) . 

[24] K. Singer et al, J. Phys. B 38, S295 (2005). 
[25] J. Stanojevic and R. Cote, Phys. Rev. A 80, 033418 
(2009). 

[26] F. Robicheaux and J.V. Hernandez, Phys. Rev. A 72, 
063403 (2005); A. Schwarzkopf, R.E. Sapiro, and G. 
Raithel, Phys. Rev. Lett. 107, 103001 (2011). 

[27] Equation Q can be derived without resorting to sym- 
metric states and operators of SA, which is however more 
involved. 

[28] A formal derivation of Eq. involves expanding the full 
derivative of g^\r) in terms of k — 0, 1, 2 photon states 
of the probe field within the blockade volume. 

[29] Our definition of the Rabi frequencies f2 PjC differ from 
that in [2l| by a factor of |. 



